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Abstract 

We consider electromagnetic field quantization in an expanding universe. We 
find that the covariant (Gupta-Bleuler) method exhibits certain difficulties when 
trying to impose the quantum Lorenz condition on cosmological scales. We thus 
explore the possibility of consistently quantizing without imposing such a con- 
dition. In this case there are three physical states, which are the two transverse 
polarizations of the massless photon and a new massless scalar mode coming 
from the temporal and longitudinal components of the electromagnetic field. An 
explicit example in de Sitter space-time shows that it is still possible to eliminate 
the negative norm state and to ensure the positivity of the energy in this the- 
ory. The new state is decoupled from the conserved electromagnetic currents, 
but is non-conformally coupled to gravity and therefore can be excited from 
vacuum fluctuations by the expanding background. The cosmological evolution 
ensures that the new state modifies Maxwell's equations in a totally negligible 
way on sub-Hubble scales. However, on cosmological scales it can give rise to a 
non-negligible energy density which could explain in a natural way the present 
phase of accelerated expansion of the universe. 

Keywords: Quantum fields in curved space-time. Dark energy 
PACS: : 04.62. -|-v, 95.36.+X 



1. Introduction 

The fact that the present phase of accelerated expansion of the universe [l| 
finds no natural explanation within known physics suggests the possibility that 
General Relativity could be inappropriate to describe gravity on cosmological 
scales. Nevertheless, despite the many attempts to find modified gravity theories 
with late time accelerated solutions 0, the models considered so far are 
generally plagued by classical or quantum instabilities, fine-tuning problems or 
local gravity inconsistencies. 

However, apart from gravity, there is another long-range interaction which 
could become relevant on cosmological scales, which is nothing but electromag- 
netism. It is generally assumed that due to the strict electric neutrality of 
the universe on large scales, the only relevant interaction in cosmology is grav- 
itation. However, this assumption seems to conflict with observations which 
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indicate that relatively strong magnetic fields with very large coherence lengths 
are present in galaxies and galaxy clusters, and whose origin is still unknown 

In this work, we consider (quantum) electromagnetic fields in an expanding 
universe. Unlike previous works, we will not follow the standard Coulomb gauge 
approach, but instead we will concentrate on the covariant quantization method. 
This method is found to exhibit certain difficulties when trying to impose the 
quantum Lorenz condition on cosmological scales (5| . In order to overcome these 
problems, and instead of introducing ghosts fields, we explore the possibility of 
quantizing without imposing any subsidiary condition. This implies that the 
electromagnetic field would contain an additional (scalar) polarization, similar 
to that generated in the massless limit of massive electrodynamics [6J] . Interest- 
ingly, the energy density of quantum fluctuations of this new electromagnetic 
state generated during inflation gets frozen on cosmological scales, giving rise 
to an effective cosmological constant 

2. Electromagnetic quantization in flat space-time 

Let us start by briefly reviewing electromagnetic quantization in Minkowski 
space-time (8|. The action of the theory reads: 

S = Jd^x (^-^F^,F^''' + A^J''^ (1) 

where is a conserved current. This action is invariant under gauge tranforma- 
tions Aft^ — >■ j4^-|-i9^A with A an arbitrary function of space-time coordinates. At 
the classical level, this action gives rise to the well-known Maxwell's equations: 

d^F^'" ^ J^. (2) 

However, when trying to quantize the theory, several problems arise related 
to the redundancy in the description due to the gauge invariance. Thus in 
particular, we find that it is not possible to construct a propagator for the 
field and also that "unphysical" polarizations of the photon field are present. 
Two different approaches are usually followed in order to avoid these difficulties. 
In the first one, which is the basis of the Coulomb gauge quantization, the gauge 
invariance of the action ^ is used to eliminate the "unphysical" degrees of 
freedom. With that purpose the (Lorenz) condition d^A^^ = is imposed by 
means of a suitable gauge transformation. Thus, the equations of motion reduce 
to: 

aA^ = Jp. (3) 

The Lorenz condition does not fix completely the gauge freedom, still it is 
possible to perform residual gauge transformations A^ — A^ -\- d^9, provided 
□0 = 0. Using this residual symmetry and taking into account the form of 
equations ([3]), it is possible to eliminate one additional component of the 
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field in tlie asymptotically free regions (typically Aq) which means ^ ■ A — 0, 
so that finally the temporal and longitudinal photons are removed and we are 
left with the two transverse polarizations of the massless free photon, which are 
the only modes (with positive energies) which are quantized in this formalism. 

The second approach is the basis of the covariant (Gupta-Bleuler) and path- 
integral formalisms. The starting point is a modification of the action in ([T]), 
namely: 

S = Jd^x (^-l^^^.F^- + |(a^A^)' + A^J^'^ . (4) 

This action is no longer invariant under general gauge transformations, but only 
under residual ones. The equations of motion obtained from this action now 
read: 

a^F^''+^a^(<9^A'') - J^. (5) 

In order to recover Maxwell's equation, the Lorenz condition must be imposed 
so that the ^ term disappears. At the classical level this can be achieved by 
means of appropriate boundary conditions on the field. Indeed, taking the four- 
divergence of the above equation, we find: 

□(a.A") = (6) 

where we have made use of current conservation. This means that the field d^A'^ 
evolves as a free scalar field, so that if it vanishes for large \t\, it will vanish 
at all times. At the quantum level, the Lorenz condition cannot be imposed 
as an operator identity, but only in the weak sense a^A-^^+^l^) = 0, where 
(-|-) denotes the positive frequency part of the operator and 10) is a physical 
state. This condition is equivalent to imposing [ao(fc) + a|| (fc)] 10) = 0, with ap 
and a|| the annihilation operators corresponding to temporal and longitudinal 
electromagnetic states. Thus, in the covariant formalism, the physical states 
contain the same number of temporal and longitudinal photons, so that their 
energy densities, having opposite signs, cancel each other. 

Thus we see that also in this case, the Lorenz condition seems to be essential 
in order to recover standard Maxwell's equations and get rid of the negative 
energy states. 



3. Covariant quantization in an expanding universe 



So far we have only considered Maxwell's theory in fiat space-time, how- 
ever when we move to a curved background, and in particular to an expanding 
universe, then consistently imposing the Lorenz condition in the covariant for- 
malism turns out to be difficult to realize. Indeed, let us consider the curved 
space-time version of action (U): 



S 



1 



-iV,,A'''f + A^J^ 



(7) 
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Now the modified Maxwell's equations read: 



and taking again the four divergence, we get: 



(8) 



(9) 



We see that once again W,yA'^ behaves as a scalar field which is decoupled from 
the conserved electromagnetic currents, but it is non-conformally coupled to 
gravity. This means that, unlike the flat space-time case, this field can be 
excited from quantum vacuum fluctuations by the expanding background in a 
completely analogous way to the inflaton fluctuations during inflation. Thus 
this poses the question of the validity of the Lorenz condition at all times. 

In order to illustrate this effect, we will present a toy example. Let us 
consider quantization in the absence of currents, in a spatially flat expanding 
background, whose metric is written in conformal time as ds^ — a{rff' (drf' — dx^) 
with a{'q) — 2-\-ta,'ah.(r] / "qt)) where 770 is constant. This metric possess two asymp- 
totically Minkowskian regions in the remote past and far future. We solve the 
coupled system of equations ([5]) for the corresponding Fourier modes of the con- 
formal field = (aAo, A), which are defined as A^{'q,x) = J d'^kA^j:{r])e'^'^^ . 

Thus, for a given mode k, the Af^ field is decomposed into temporal, longitudinal 
and transverse components. The corresponding equations read: 



A" 



Ofe 



J - 2H' 



Aah - 2ik 



1 + e 



2? 



-A'a 



■Ok 



2e 

HAok 



A'„k - HA 



life 



= 



k^Ai 







(10) 



withH — a'/aandfc ~ \k\. We see that the transverse modes are decoupled from 
the background, whereas the temporal and longitudinal ones are non-trivially 
coupled to each other and to gravity. Let us prepare our system in an initial 



state belonging to the physical Hilbert space, i.e. satisfying d^A. 







in the initial flat region. Because of the expansion of the universe, the positive 
frequency modes in the in region with a given temporal or longitudinal polar- 
ization A will become a linear superposition of positive and negative frequency 
modes in the out region and with different polarizations A' (we will work in the 
Feynman gauge ^ = —1), thus we have: 



E 

A'=0, 



[a,A'(fc)^'V°"*^ 



fi k 



l3xx'{k)A 



X' (out) 
^ —k 



or in terms of creation and annihilation operators: 
al™*)(fc) 



E 

A=0, 



[a,v(fc)afHfc) 



(11) 



(12) 
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with A, A' = 0,11 and where a\y and /3aa' are the so-called Bogoliubov co- 
efficients (see [9| for a detailed discussion), which are normalized in our case 
according to: 

("Ap "A'p' Vpp' - PxpPx'p' Vpp' ) = Vw (13) 

P.P'=0,|| 

with rjw = diag{— 1,1) with A, A' — 0, ||. Notice that the normalization is 
different from the standard one 9], because of the presence of negative norm 
states. 




Figure 1: Occupation numbers for temporal (continuous line) and longitudinal (dashed line) 
photons in the out region vs. k in rj^^ units. We also show the analytical approximate solution 
(dotted line) given in the main text. 

Thus, the system will end up in a final state which no longer satisfies the 
weak Lorenz condition i.e. in the out region d^Ag^^''\(j)) ^ 0. This is shown in 
Fig. 1, where we have computed the final number of temporal and longitudinal 
photons n'y*{k) = l/?AA'(fc)P, starting from an initial vacuum state with 
ng"(A;) = n||"(fc) = 0. We see that, as commented above, in the final region 
nQ"*(fc) 7^ n||"*(fc) and the state no longer satisfies the Lorenz condition. Notice 
that the failure comes essentially from large scales {kijo 1), since on small 
scales {krjo ^ 1), the Lorenz condition can be restored. This can be easily 
interpreted from the fact that on small scales the geometry can be considered 
as essentially Minkowskian. 

We can obtain an approximate analytical solution for the Bogoliubov coeffi- 
cients by assuming an expansion in which H vanishes in the in and out regions, 
but it takes a constant value h in the period {—to, to). Moreover, we shall ne- 
glect the production of longitudinal modes (which is justified from our numerical 
computation). For this simplified problem, the relevant Bogoliubov coefficients 
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are given by: 



Bun - ^2to(h+ik) 



In'' 



hi 

h 

h J ■ y/Ah^ ~ fc2 



— cosh uj H , smh t 



ik 



2aLt h 



2/i( 



l2 „V4/.2 _ P 



sinhoj (14) 



where w = 2tQ\/ AK^ — fc^, ao„f is the scale factor in the out region and we 
have assumed that ai„ = 1. In Fig. 1 we show the corresponding production 
of temporal modes = |/3ooP + |/3||oP for h = 0.29 and io = 2 (which 

approximately corresponds to the Hubble expansion rate used for the numeri- 
cal computation). Notice that the approximate analytical solution reproduces 
very accurately the numerical result at small k. This is because, on that re- 
gion, the production of longitudinal modes is practically zero and for very large 
wavelenghts, the production of temporal states is not sensitive to the particular 
shape of the expansion rate. 

These results show that in a space-time configuration with asymptotic flat 
regions, an initial state satisfying the weak Lorenz condition does not necessarily 
satisfy it at a later time. In order to overcome this problem, it is possible to 
impose a more stringent gauge-fixing condition. Thus, following [lO| we can 
define the physical states \(j)) as those such that V^y4.^*^+^ |(/)) = 0, Vr;. Although 
this is a perfectly consistent solution, notice that the separation in positive and 
negative frequency parts depends on the space-time geometry and therefore, 
the determination of the physical states requires a previous knowledge of the 
geometry of the universe at all times. 

Another possible way out would be to modify the standard Gupta-Bleuler 
formalism by including ghosts fields as done in non-abelian gauge theories [ll| . 
With that purpose, the action of the theory ^ can be modified by including 



the ghost term (see [12|): 

Sg^ [ d^x^g^''d,,cd,c (15) 



where c are the complex scalar ghost fields. It is a well-known result 
that by choosing appropriate boundary conditions for the electromagnetic and 
ghosts Green's functions, it is possible to get {(j)\T^^ + Tl^^\(f>) — 0, where T^^ 
and denote the contribution to the energy-momentum tensor from the ^ 
term in ([7]) and from the ghost term p5)l respectively. Notice that a choice of 
boundary conditions in the Green's functions corresponds to a choice of vacuum 
state. Therefore, also in this case an a priori knowledge of the future behaviour 
of the universe geometry is required in order to determine the physical states. 

In this work we follow a different approach in order to deal with the diffi- 
culties found in the Gupta-Bleuler formalism and we will explore the possibility 
of quantizing electromagnetism in an expanding universe without imposing any 
subsidiary condition. 
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4. Quantization without the Lorenz condition 



In the previous section it has been shown that ahhough the Lorenz gauge- 
fixing conditions can be imposed in the covariant formalism, this cannot be 
done in a straightforward way. These difficulties could be suggesting some more 
fundamental obstacle in the formulation of an electromagnetic gauge invariant 
theory in an expanding universe. As a matter of fact, electromagnetic models 
which break gauge invariance on cosmological scales have been widely considered 
in the context of generation of primordial magnetic fields (see, for instance, [13] )• 

Let us then explore the possibility that the fundamental theory of electro- 
magnetism is not given by the gauge invariant action ([1]), but by the gauge 
non- invariant action Notice that although ([7]) is not invariant under gen- 
eral gauge transformations, it respects the invariance under residual ones and, 
as shown below, in Minkowski space-time, the theory is completely equivalent 
to standard QED. Since the fundamental electromagnetic theory is assumed 
non-invariant under arbitrary gauge transformations, then there is no need to 
impose the Lorenz constraint in the quantization procedure. Therefore, having 
removed one constraint, the theory contains one additional degree of freedom. 
Thus, the general solution for the modified equations (H)) can be written as: 

A^ = Ai'^+A^^^+A^^^+d^9 (16) 

where A^/l"^ with i = 1,2 are the two transverse modes of the massless photon, 
A'lf ^ is the new scalar state, which is the mode that would have been elim- 
inated if we had imposed the Lorenz condition and, finally, df^O is a purely 
residual gauge mode, which can be eliminated by means of a residual gauge 
transformation in the asymptotically free regions, in a completely analogous 
way to the elimination of the Aq component in the Coulomb quantization. The 
fact that Maxwell's electromagnetism could contain an additional scalar mode 
decoupled from electromagnetic currents, but with non-vanishing gravitational 
interactions, was already noticed in a different context in 

In order to quantize the free theory, we perform the mode expansion of the 
field with the corresponding creation and annihilation operators for the three 
physical states: 

A,^fd'k J2 hA(fc)^il^+at(fc)^] (17) 

where the modes are re quir ed to be orthonormal with respect to the scalar 
product (see for instance |10|): 



i I dS^ 

s 



= SxyS'^'^Hk-k'), A,A' = l,2,s (18) 



where dS^ is the three-volume element of the Cauchy hypersurfaces. In a 
Robertson- Walker metric in conformal time, it reads dS^ = a'^{ri){cfix, 0,0,0). 
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The generalized conjugate momenta are defined as: 



(19) 



Notice that the three modes can be chosen to have positive normahzation. The 
equal-time commutation relations: 



[A^i7j,x),AAv,x')] = [nO'^(7?,f),nO'^(r;,f/)] =0 



and 



(20) 



(21) 



(22) 



by means of the normalization condition in (|T8|) . Notice that the sign of the 
commutators is positive for the three physical states, i.e. there are no negative 
norm states in the theory, which in turn implies that there are no negative 
energy states as we will see below in an explicit example. 

Since V fj_A^ evolves as a minimally coupled scalar field, as shown in ([S]), on 
sub-Hubble scales (|fc?7| ^ 1), we find that for arbitrary background evolution. 



[^^(77,x),n<'^(r?,f' )] = z^S^'H^-x') 
can be seen to imply the canonical commutation relations: 

aA(fc),al,(fc')l =SxyS(^^\k^k'), A,A' = 1,2,. 



W^A 



(s)M| 



cx a 



i.e. the field is suppressed by the universe expansion, thus 



effectively recovering the Lorenz condition on small scales. Notice that this is 
a consequence of the cosmological evolution, not being imposed as a boundary 
condition as in the flat space-time case. 

On the other hand, on super-Hubble scales {\kr]\ ^ 1), |V^^^^''^| = const. 
which, as shown in Q, implies that the field contributes as a cosmological con- 
stant in (O. Indeed, the energy-momentum tensor derived from ([7]) reads: 



T — 



(Vo^")' -t- 2A"V„ (VpA^)] - 4A(^V,) (VoA")l (23) 



Notice that for the scalar electromagnetic mode in the super-Hubble limit, the 
contributions involving F^;^ vanish and only the piece proportional to ^ is rele- 
vant. Thus, it can be easily seen that, since in this case V^A" = constant, the 
energy-momentum tensor is just given by: 



.9p.(V„A")2 



(24) 



which is the energy- momentum tensor of a cosmological constant. Notice that, 
as seen in the new scalar mode is a massless free field. This is one of the 
most relevant aspects of the present model in which, unlike existing dark energy 



8 



theories based on scalar fields, dark energy can be generated without including 
any potential term or dimensional constant. 

Since, as shown above, the field amplitude remains frozen on super-Hubble 
scales and starts decaying once the mode enters the horizon in the radiation or 
matter eras, the effect of the ^ term in ([5]) is completely negligible on sub-Hubble 
scales, since the initial amplitude generated during inflation is very small as we 
will show below. Thus, below 1.3 AU, which is the largest distance scale at 
which electromagnetism has been tested [l^, the modified Maxwell's equations 
([H]) are physically indistinguishable from the flat space-time ones ^ . 

Notice that in Minkowski space-time, the theory ([T]) is completely equivalent 
to standard QED. This is so because, although non-gauge invariant, the corre- 
sponding effective action is equivalent to the standard BRS invariant effective 
action of QED. Thus, the effective action for QED obtained from ([2]) by the 
standard gauge-fixing procedure reads: 

e'^ = J [d^][dc][dc][#][#]e'^'''^(-^^-^''''+^(^'*^'')'+'''''^''^"^''^+^^) (25) 

where Cp is the Lagrangian density of charged fermions. The ^ term and the 
ghosts field appear in the Faddeev-Popov procedure when selecting an element 
of each gauge orbit. However, ghosts being decoupled from the electromagnetic 
currents can be integrated out in fiat space-time, so that up to an irrelevant 
normalization constant we find: 

e^^ = I [dA][dV'][#]e'/'^'-(^^^-^"''+^(^-^^)'+^-) (26) 

which is nothing but the effective action coming from the gauge non-invariant 
theory ((T]) in fiat space-time, in which no gauge-fixing procedure is required. 



5. An explicit example: quantization in de Sitter space-time 



Let us consider an explicit example which is given by the quantization in 
an infiationary de Sitter space-time with 0(77) = — l/(_ff/7y). Here we will take 
^ — 1/3, although similar results can be obtained for any f > 0. For this case, 
in which H' = H^, we can obtain the following expression for the temporal 
component in terms of the longitudinal one from (|10|) : 



A, 



Ok 



2fc(4fc2 + 3-H2 



d^A 



life 



drj^ 



n- 



d^A 



life 



drj'^ 



drj 



life 



(27) 



On the other hand, if we insert this relation into the equation for A\\k we obtain 
a fourth-order equation for the longitudinal component which is completely 
decoupled from A^k- This fourth-order equation can be expressed as follows: 



.2P 



d 

3^2 



4fc4 + - 9k^H^ 
4fc2 + 3-H2 



d^ 

diff 



2H 



drj 



A\\k = 0, 
(28) 
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where we have performed a factorization in such a way that the second bracket 
is nothing but the operator determining the evolution of a free scalar field, i.e., 
D9 = 0. Now, we shall use the notation J^{Q{A\\k)) — 0, where J" and G are 
the operators given in the first and second brackets in (^5)) respectively. In 
particular, ^ = □ as we have already said. Thus, in order to solve the equation 
we have to obtain the kernel of J", whose solutions will be denoted by S, and, 
then, solve the equation for a free scalar field with S as an external source: 
G{A\\k) = S. This method has the advantage of allowing to identify the pure 
residual gauge mode as the homogeneous solution of this last equation. The 
explicit solution for the normalized scalar state is: 



1 



(27r)3/2 y2fc 



1 

krj 



-(1 + ikr])e~'^'^ ~ k^ri^e''''^Ei{2ikr)) 



Akx 



A 



(^) 
life 



{(l + z/c77)e-''=" 



.e-»fer, ^ _ ih])e''''^ Ei{2ikT]) 



Akx 



(29) 



where Ei{x) = e~'*^ /tdt is the exponential integral function and we have 
fixed the integration constants so that the mode is canonically normalized ac- 
cording to (fT8| . Using this solution, we find: 



\vi) ik 3 (1 + ifcr?) 



kr]-\-ikx 



(27r)3/2V2fc2 A:2r;2 



(30) 



so that the field is suppressed in the sub-Hubble limit as V^^^,*'*'^ ~ 0{{kr])~'^). 

On the other hand, from the energy density given by = T^q, we obtain 
in the sub-Hubble limit the corresponding Hamiltonian, which is given by: 



k J2 [al(fc)a;,(fc)+a;,(fc)aUfc) 



\=1.2,s 



(31) 



We see that the theory does not contain negative energy states (ghosts). In fact, 
as shown in the theory does not exhibit either local gravity inconsistencies 

or classical instabilities. 

Finally, from ((29)) it is possible to obtain the dispersion of the effective 
cosmological constant during inflation: 



rik 

{0\iV,An''\0) = I -jPAik) 



(32) 



with PA{k) = 47rfc3|V,,yl[^^^^|2. In the super-Hubble limit, we obtain for the 
power-spectrum: 



PA{k) = 



9Hj 
167r2 



(33) 
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in agreement with Notice that this result imphes that pA ~ (Hj)'^. The 
measured value of the cosmological constant then requires Hj ~ 10"'^ eV, which 
corresponds to an inflationary scale Mj ^ 1 TeV. Thus we see that the cosmo- 
logical constant scale can be naturally explained in terms of physics at the 
electroweak scale. 

Once we have computed the initial amplitude of the electromagnetic fluctu- 
ations, we can estimate the magnitude of the corrections induced by the new 
terms in the modified Maxwell's equations. Notice that ([SJ can be rewritten as 
the ordinary Maxwell's equation with an additional conserved current source: 

W^F'"' = Ji' + (34) 

where = -^S/f^i^/ ^A") satisfies V^J^^ = by virtue of dH). For sub-Hubble 
modes below 1.3 A.U., which is the region on which electromagnetism has been 
tested, consistency requires the new current to be negligible. It is possible to 
estimate the present value of such a current by noticing that a Fourier mode k 
entered the Hubble radius in the radiation dominated era when the scale factor 

1 /2 

was ain — Hq aeq /fc, with Ueq the scale factor at matter-radiation equality and 
where we have taken that today oq = 1. Therefore, since the initial amplitude 
V^A'' ~ remains constant until horizon reentry, today we obtain for the 
effective electric charge density created by the mentioned modes: 

J| - fc iJ| a^n - Ho Hj a)l^ - lO"'^^ eV^ (35) 

which is independent of k and completely negligible since it roughly corresponds 
to the charge density of one electron in the volume of the Earth. 

6. Conclusions 

In this work we have considered the difficulties in the application of the 
covariant quantization method for electromagnetic fields in an expanding uni- 
verse. Although the Lorenz gauge-fixing condition can in principle be formally 
applied also in cosmological contexts, we have explored the alternative possi- 
bility of quantization without subsidiary conditions, and therefore, without the 
introduction of ghost fields. 

In this approach the theory is invariant only under residual gauge transfor- 
mations and a new scalar mode appears for the electromagnetic field. Despite 
these facts, standard QED is recovered in the flat space-time limit as the new 
scalar state is completely decoupled from the conserved electromagnetic currents 
(only transverse photons couple to them [l3l ) . The residual gauge symmetry of 
the theory ensures that the negative norm state can be eliminated and that the 
energy is positive as shown in an explicit example. 

This quantization procedure is found to have interesting consequences for 
dark energy. Thus, the energy density of the new mode on super-Hubble scales 
(which is essentially the temporal part of the electromagnetic field) is seen to 
behave as a cosmological constant irrespectively of the expansion rate of the 
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universe. The new mode can be generated during inflation, in a similar way 
to the inflaton fluctuations and this fact aUows to establish a link between the 
scale of inflation and the value of the cosmological constant. 
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